Inclusive nucleon emission induced by quasi-elastic neutrino-nucleus interactions 
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A previous model on inclusive charged-current quasi-elastic nuclear reactions Q is extended to 
include neutral- and charged-current nucleon emission reactions. The problem of outgoing nucleon 
propagation is treated by means of a Monte Carlo simulation. 



I. INTRODUCTION 

The main effort of theoretical studies on neutrino- 
nucleus reactions at intermediate energies (few hundreds 
MeV) has been devoted to quasi-elastic (QE) charged- 
current (CC) reactions. This is due to the fact that water 
Cerenkov detectors only detect outgoing charged leptons. 
New experiments such as K2K aim to measure the espec- 
trum of outgoing nucleons produced in neutral currents 
(NC) in this energy range. These kind of observables are 
also of interest when dealing with the extraction of in- 
formation on strange quark axial content of the nucleon. 
For a correct theoretical anlysis of these experiments we 
have to take into account the rescattering effect of the 
outgoing nucleons in its way out of the nucleus. 

In [l| a model dealing with QE CC reactions of the type 
J ^X(vi,l~)Y was presented. This model was based in a 
local Fermi gas model of the nucleus and included further 
nuclear effects in a parameter free fashion, thus leading to 
one of the best descriptions of the LSND experiments. 
Our confidence in this model at intermediate energies 
higher than those of is based on the fact that a similar 
model [3] is the only model on electron scattering that 
has been successfully compared with experiment from the 
QE to the A excitation regions. 

In section [n] we extend the model in [jj to NC 
zX(v h i>i)Y . Section |Inl deals with the Monte Carlo 
treatment of the rescattering problem. Finally we show 
results for outgoing nucleon observables in CC and NC. 



II. NEUTRAL CURRENTS 

In NC processes the Z° boson can be absorbed by one 
nucleon leading to the QE contribution of the nuclear 
response function. 
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We obtain the cross section of this proccess by calculating 
the imaginary part of the self-energy diagram associated 
to the Z° absorption in nuclear matter. Results for finite 
nuclei are obtained with the local density approximation. 



This method was also used in pj with CC processes, the 
main differences in NC arising from the vertex Z°NN 
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with vector and axial nucleon NC given by 
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Thanks to SU(3) symmetry some relations exist among 
the NC form factors and the CC (F^ 2 , G A and G P ) 
and the electromagnetic ones. We have to introduce 
F" , fi s Fi ,G S A and G S P , the strange vector and axial nu- 
cleon form factors Q| . We use the parametrizations of [3 
for electromagnetic form factors and those of fit II in [6j 
for strange ones. 

Like in [l| further nuclear effects were taken into ac- 
count, mainly keeping a correct energy balance and in- 
cluding long range correlations (RPA). The main dif- 
ferences with respect to CC arise from the inclusion of 
isoscalar pieces in the nucleon-nucleon effective interac- 
tion and the absence of the outgoing lepton Coulomb 
interaction with the final nucleus. 



III. MONTE CARLO SIMULATION 

In the previous section we overviewed the evaluation of 
NC inclusive neutrino induced cross sections in the QE 
region. Thus we determined, for a fix incoming neutrino 
energy, the inclusive QE cross section cPa/dQ'dE' (CI', 
E' solid angle and energy of the outgoing lepton). In 
our scheme this is obtained after performing an integra- 
tion over the whole nuclear volume, following the line of 
the local density approximation of the Fermi gas model. 
Thus, for a fix transferred four momentum chosen ac- 
cording to d 2 a I dQ! dE' , we can randomly select the point 
of the nucleus where the absorption takes place using the 
profile d 5 cr/dn'dE'd 3 r. 

Now, we need to have on top of that the distribution 
of three-momenta of the outgoing nucleon. Since events 
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are generated probabilistically one by one, we generate 
a random momentum from the local Fermi sea, and this 
will generate an outgoing momentum p' = p+q. If it hap- 
pens that |p'| < kp(r) (the local Fermi momentum) then 
the event is Pauli blocked, it is dismissed and another 
event is generated. Thus, we have already the configu- 
ration of the final state after the first step: in this case 
just one nucleon produced in the point r of the nucleus 
with momentum p'. As with respect to having a proton 
and a neutron in the final state, this is trivially done: for 
the CC case, the outgoing nucleon is a proton (neutron) 
for neutrino (antineutrino) induced processes, while for 
NC the reaction probability was already splitted into a 
proton and a neutron induced ones. 

Then we simulate the trajectory of the ejected nucleon 
in its way out of the nucleus. We move the nucleon 
in finite length steps (d << A, the nucleon mean free 
path) . At every new position we evaluate the probability 
of the nucleon colliding so it produces a new outgoing sec- 
ondary nucleon. All of the ejected nucleons (the primary 
one, and the possibly existing secondary) are again prop- 
agated until they all leave the nucleus, all of them giving 
the same contribution to the differential cross section. 

When a QE nucleon-nucleon collision is produced a 
random Fermi sea nucleon momentum is selected. The 
new primary and secondary nucleon kinematics are se- 
lected according to experimental differential cross sec- 
tions of free space nucleon-nucleon collision, conveniently 
modified in order to take into account the effects of 
medium polarization (correlations) and Pauli blocking 
effects. The propagation of nucleons is simulated with 
a semiclassical approximation of nucleons moving in 
a nucleon-nucleus potential. This approximation was 
shown to be fairly accurate for nucleons with kinetic en- 
ergies higher than ~ 20 MeV in 0. 




FIG. 2: Neutral current 40 O(u, u + N) at 500 MeV (upper panels) 
and 150 MeV (lower panels) cross sections as a function of the 
kinetic energy of the final nucleon. Left and right panels correspond 
to the emission of protons and neutrons respectively. The solid 
histogram shows results without FSI and the dashed one the full 
model. 



form factor of the nucleon, and thus to the g s A parame- 
ter. Our rescattering approach produces minor changes 
for light nuclei, because of the smaller average density, 
and for low energies because most secondary nucleons 
are below our 30MeV cut. However, the sensitivity is 
larger for both heavier nuclei and for larger energies of 
the neutrinos as shown in Fig. |3 



IV. RESULTS 

The nucleons spectra produced by CC processes in- 
duced by muon neutrinos are shown in Fig. ^ for Argon. 
Of course neutrinos only interact via CC with neutrons 
and would emit protons, but these primary protons in- 
teract strongly with the medium and collide with other 
nucleons which are also ejected. As a consequence there 
is a reduction of the flux of high energy protons but a 
large number of secondary nucleons, many of them neu- 
trons, of lower energies appear. 

The energy distributions of nucleons emitted after a 
NC interaction are shown in Fig. [3 Our results with- 
out rescattering compare well with other calculations like 
those of [U- However, in these latter cases, which incor- 
porate the nucleons final state interaction via the use of 
optical potentials, the main effect is to reduce the cross 
section at all energies instead of displacing the strength 
towards lower energies as we find. 

The ratio of proton to neutron quasielastic cross sec- 
tion could be very sensitive to the strange quark axial 
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FIG. 3: Ratio of da/dE for protons over that for neutrons for 
E y = lhOMeV and E v = 500MeV in the reaction v + 40 At -> 
v 1 + N + X as a function of the nucleon kinetic energy. 
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